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Based on analyses of multiple TCR transgenic (tg) models, the emergence of pathogenic T cells 
in diabetes-prone NOD mice has been ascribed to a failure to censure autoreactive clones in the 
thymus. In contrast, using isolated and preselected thymocytes, we show that nonobese diabetic 
(NOD) genetic variation impairs neither clonal deletion nor downstream transcriptional programs. 
However, we find that NOD genetic variation influences ap/78-lineage decisions promoted by 
early expression of tg ap-TCRs at the double-negative (DN) stage. In BG and other genetic 
backgrounds, tg ap-TCRs behave like 78-TCRs and commit a large fraction of DNs toward 
the 78-lineage, thereby decreasing the size of the double-positive (DP) pool, which is efficiently 
positively and negatively selected. In NOD DNs, ap-TCR signalosomes instead behave like pre- 
TCRs, resulting in high numbers of DPs competing for limited selection niches, and poor positive 
and negative selection. Once niche effects are neutralized in mixed bone marrow chimeras, 
positive and negative selection are equally efficient on B6 and NOD backgrounds. Biochemical 
analysis revealed a selective defect in the activation of Erk1/2 downstream of NOD aP-TCR 
signalosomes. Therefore, NOD genetic variation influences ap/78-lineage decisions when the 
ap-TCR heterodimer is prematurely expressed, but not the process of negative selection. 



Two major mechanisms contribute to the estab- 
lishment of a self-restricted but mostly non- 
autoreactive T cell repertoire (Starr et al., 2003). 
First, after the onset of a(3-TCR expression in 
CD4 + CD8 + (double positive [DP]) thymocytes, 
positive selection rescues from programmed cell 
death those cells expressing TCRs that recognize 
with sufficient (but low) affinity self-peptides 
complexed to MHC molecules. As the second 
arm of central tolerance, negative selection — 
mainly through clonal deletion — eliminates from 
the repertoire T cells that are equipped with 
TCRs of higher affinities (Siggs et al., 2006), 
although a small subset of cells expressing self- 
reactive TCRs is rescued and diverted toward 
regulatory lineages such as regulatory T cells 
(T reg cells; Baldwin et al., 2004). 

Mechanistically, the opposite biological out- 
comes — survival and cell death — promoted by 
the same a(3-TCR in conditions of positive and 
negative selection are driven by temporally, spa- 
tially, quantitatively, and qualitatively distinct 
signals. One of the most differential pathways is 
the RasGRPl -driven activation ofErkl/2 MAP 
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kinases (Teixeiro and Daniels, 2010), which plays 
an obligatory role during positive selection but is 
dispensable for negative selection (McGargill 
et al., 2009). Conversely, Grb2-SOS— driven ac- 
tivation of p38 and JNK MAP kinases is required 
for negative but not positive selection (Siggs 
et al., 2006). A few other molecules are known to 
participate in negative selection, including Erk5, 
Bim, Nur77, and Mink (Sohn et al., 2007). 

Although there is ample evidence for the 
importance of dominant tolerance andT reg cells 
for preventing autoimmunity, data supporting an 
equivalent role for deletional tolerance are less 
abundant. Some evidence has come from the auto- 
immune disease caused by the disruption of the 
Aire gene (Mathis and Benoist, 2009), and a role 
for central tolerance defects in the pathogenesis 
of type-1 diabetes (T1D) has also been proposed, 
based on several observations. First, genetically 
controlled variation in the level of insulin gene 
expression in the human thymus has been linked 
to susceptibility to T1D, and proposed to reflect 
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a variable degree of tolerance induction (Pugliese et al., 1997). 
Second, nonobese diabetic (NOD) thymocytes are resistant to 
anti-CD3— provoked clonal deletion, although this point has 
been debated (Kishimoto and Sprent, 2001; Villunger et al., 
2003). Several TCR transgenic (tg) derivatives of the NOD 
mouse model of T1D have been used to address this issue. 
In particular, studies using the diabetogenic BDC2.5 and AI4 
TCRs that recognize natural autoantigens in the pancreas 
(Choisy-Rossi et al, 2004; Zucchelli et al., 2005; Holler et al, 
2007; Serreze et al, 2008) or the 3A9 TCR specific for a 
transgene-encoded pancreatic neo-antigen (Lesage et al., 2002; 
Liston et al., 2004, 2007) have shown that autoreactive clones 
appear not to be purged efficiently in NOD mice. Several of 
the faulty molecular candidates were found to be common to 
all three TCR tg models, including Bim and Nnr77 genes. 
A body of literature, not easily integrated, indicates that these 
defects may result from proximal and multifocal signaling 
defects downstream from the TCR in NOD mice (Rapoport 
et al, 1993, 1999; Salojin et al, 1997; Zhang et al, 1998). 

A challenge remains to understand how signaling and 
transcriptional defects downstream of a(3-TCRs on the NOD 
background lead to selective functional outcomes, and how 
these are affected by NOD genetic variation. In the models 
used previously by ourselves and others, careful dissection was 
potentially confounded by effects of NOD genetic variation 
at multiple stages of thymocyte differentiation. Here, by using 
experimental systems in which this parameter was rigorously 
controlled, we arrived at a different conclusion: the phenom- 
ena interpreted as resistance to clonal deletion on the NOD 
background result instead from increased competition for se- 
lection niches, itself a consequence of a differential ability 
of NOD cells to proceed into ct(3 over 78 commitment in 
response to a(3-TCR expression at the immature double- 
negative (DN) stage. This early divergence in TCR tg models 
reveals an unanticipated impact of NOD genetic variation on 
a crucial node in a(3-TCR signaling networks, Erkl/2. 

RESULTS 

Clonal deletion and associated transcriptional programs 
are functional In NOD mice 

Our earlier results (Zucchelli et al., 2005; Holler et al., 2007) 
showed that, as with other TCR tg systems (Lesage et al., 
2002; Liston et al., 2004; Choisy-Rossi et al, 2004; Serreze 
et al, 2008), expression of the BDC2.5 TCR led to a dearth 
of DPs and small thymi on the B6 g7 but not on the NOD 
background, observations which were interpreted as a differ- 
ential resistance to clonal deletion. Thymocyte differentiation 
is a dynamic process involving several TCR-dependent check- 
points on which genetic variations can impinge. We reasoned 
that monoclonal and synchronized DP thymocytes, stimulated 
in vitro with their cognate ligand, would allow us to identify 
how genetic variation has an impact on the TCR signaling 
axis of DP thymocytes while strictly controlling other vari- 
ables (in particular, potential effects of NOD genetic variation 
at other stages of thymocyte differentiation). Thus, we bred 
the BDC2.5 TCR transgene onto the B6.Rag _/ " and NOD. 
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H2 b .Rag~ 7 ~ backgrounds, both expressing the nonselecting 
A b MHC molecule. As anticipated, the resulting mice did not 
harbor mature T cells and thymocyte differentiation was 
arrested at a preselection DP stage (Fig. SI, A and B). 

Purified DP thymocytes from these mice were stimulated 
in vitro with mimotope-pulsed splenic APCs, and early conse- 
quences on the transcriptome were analyzed by gene expression 
profiling. A robust stimulation was observed after 7 h. Contrary 
to anticipated results, however, the responses were very similar 
on the NOD.H2 b and B6 backgrounds, with most of the re- 
sponsive genes aligning on the diagonal of the fold change/fold 
change plot (Fig. 1 A), including at shorter time points and with 
suboptimal peptide concentrations (Fig. SI C). This finding 
applied not only to T cell activation genes (Fig. 1 B) but also to 
negative selection— associated (Nr4al, Bcl2ll 1, Gadd45b, and Pdcd'l) 
and prosurvival (Rorc, Hivep2, Bcl2U, and Bel 2) genes (Fig. 1 C), 
many of which were previously found to be improperly regu- 
lated during negative selection in NOD mice (Choisy-Rossi 
et al, 2004; Liston et al, 2004, 2007; Zucchelli et al, 2005). 

We also asked whether varying the unit affinity of the li- 
gand might reveal a difference, using mimotope variants with 
different agonist potencies (Fig. 1 D). CD69 up-regulation, 
TCR down-modulation, and ligand-induced cell death were as 
efficient in NOD.H2 b as in B6 thymocytes, through a range of 
mimotope doses and potencies. In conclusion, preselected DP 
thymocytes on the NOD background showed largely unaltered 
TCR-dependent transcriptional and apoptotic responses. 

Positive selection of class I- and class II— restricted tg TCRs 
is impaired in NOD mice 

These results disagree with earlier conclusions made by ourselves 
and others, reporting central tolerance defects in NOD thymi 
in vivo or in fetal thymic organ cultures (FTOCs; Lesage et al., 
2002; Choisy-Rossi et al, 2004; Liston et al, 2004; Zucchelli 
et al., 2005). In these systems, the NOD background was 
shown to impede the deletion of autoreactive DP and/ or CD4 
single-positive (SP) thymocytes, resulting in accumulation of the 
corresponding populations. Because whole thymi, rather than 
isolated cell populations, were exposed to negatively selecting li- 
gands in these studies, compound effects of NOD alleles at sev- 
eral stages of thymocyte differentiation might have confounded 
the results. We thus decided to revisit the impact of NOD 
genetic variation on thymocyte selection, using adult BDC2.5 
TCR tg mice on selecting NOD or B6 g7 backgrounds. 

Similar to what was observed at the fetal stage (Zucchelli 
et al., 2005), thymi from adult .Rag-sufficient and -deficient 
BDC2.5/B6 g7 mice showed a markedly reduced proportion of 
DPs, a much higher proportion of DNs, more CD4SPs, and 
lower total cellularity (Fig. 2, A and B) relative to BDC2.5/NOD 
thymi, whose overall appearance and cellularity were comparable 
with those of a non-tg thymus. This dearth of DPs relative to SPs 
and DNs could be interpreted as reflecting various degrees of 
negative selection, yet they were also compatible with a differ- 
ence in efficacy of positive selection and of the DP to SP transi- 
tion. Indeed, the following careful analysis of these thymi 
supported a defect in positive selection in BDC2.5/NOD mice: 
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Figure 1. Negative selection is functional in NOD mice. DP thymocytes purified from mice expressing the BDC2.5 tg TCR on Rogl-deficient, 
nonselecting B6 or N0D.H2 b backgrounds, were co-cultured with mimotope-pulsed splenic APCs. DP thymocytes were analyzed by transcriptome profiling 
(A-C) or flow cytometry (D). (A) Fold change/fold change plot showing transcriptome changes induced by a 7-h stimulation with 500 nM BDC mimotope 
peptide, as the stimulated (stim.J/unstimulated (unstim.) ratio of expression in B6 (x axis) and NOD (y axis) DP thymocytes. For the unstimulated condition, 
DP thymocytes were incubated with unpulsed splenic cells. (B) Heat map of the 50 transcripts most induced in DPs by mimotope stimulation, in the same 
conditions as in A. Relative expression is color coded for two peptide concentrations and three time points (left margin). For the unstimulated (unstim.) 
condition, DP thymocytes were incubated with unpulsed splenic cells for 3 h. (C) Microarray expression values for selected proapoptotic and prosurvival 
genes in unstimulated or stimulated cells as in A, plotted versus time. (D) Stimulation by three mimotope variants (BDCmi#1-3) of different potency. 
Frequency of CD69+ cells among DP thymocytes (top), geometric mean of CD3e expression (middle), and ligand-induced apoptosis, measured as the fre- 
quency of viable DPs (bottom), after a 20-h stimulation. Error bars indicate standard deviation of the mean. Data are means of two (A-C; n = 2 per genotype 
for t = 1 and 3 h and n = 4 per genotype for t = 7 h) and representative of three (D; n > 2 per genotype per experiment) independent experiments. 



(A) The frequency of CD4SP thymocytes and peripheral 
T cells with a lower expression of the BDC2.5 clonotypic TCR 
(BDC low population) was three to four times higher on the 
NOD background (Fig. 2, C and D).This resulted from increased 
rearrangement of endogenous genes encoding Vet chains and 



surface expression of two TCR-a chains on the cell surface be- 
cause BDC low CD4 + cells were absent in R^-deficient BDC2.5 
mice (Fig. S2). Endogenous rearrangements of the Tern locus were 
reported to be promoted in TCR tg models by defects in positive 
selection, rather than negative selection (Buch et al, 2002). 
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Figure 2. Impaired positive selection of BDC + thymocytes in NOD mice. Thymocytes were isolated from polyclonal or BDC2.5 TCR tg mice on B69 7 or NOD 
backgrounds. (A) Distribution of thymocyte subsets. (B) Number of total thymocytes, versus the age of each mouse. (C) Expression of the BDC2.5TCR, detected 
with the anti-clonotypic BDC antibody (Kanagawa et al, 2002), on CD4SP thymocytes (top) and CD4 + splenocytes (bottom). (D) Quantification of data shown in C. 
Horizontal lines indicate the mean. (E and F) Expression of CD3e on DP and CD4SP thymocytes from BDC2.5TCR tg (E) and polyclonal mice (F). Thin lines indicate 
negative control staining. (G) Percentage of BDC low CD4SP thymocytes plotted versus the geometric mean of CD3e surface expression on DP thymocytes from mice 
of indicated genotypes. (H) BrdU pulse-labeling experiment. Number of BrdU + CD4SP per 1000 BrdU + DP thymocytes over time, after a single BrdU injection. Error 
bars indicate standard deviation of the mean. Data are from 15 (A-G; n > 40 per genotype) and means of 8 (H; n> 25 per genotype) independent experiments. 



(B) DP thymocytes from BDC2.5/NOD mice expressed 
40-50% less CD3 than did their B6e 7 counterparts (Fig. 2 E). 
This difference was absent in polyclonal mice (Fig. 2 F) and, 
thus, did not reflect a general defect in TCR expression. TCR 
levels on DP thymocytes from TCR tg models positively cor- 
relate with the efficiency of positive selection (Sato et al., 
1996), a regulation which applied to the BDC2.5 mouse model 
because the gradual decrease of the availability of selecting A B? 
molecules led to a commensurate decrease of TCR levels on 
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the B6 background (Fig. 2 G). On this plot, BDC2.5/NOD 
DPs were roughly analogous to those inefficiently selected on 
the hybrid BDC2.5/B6s 7/b background. 

(C) We directly evaluated the efficacy of positive selection 
by performing BrdU pulse-labeling experiments. This tracer is 
incorporated by cycling DN and early DP thymocytes and 
only appears in CD4SPs a few days later, with their rate of ap- 
pearance reflecting the efficacy of positive selection. Equiva- 
lent proportions of thymocytes were labeled during the pulse 
period (Fig. S3), but the emergence of BrdU + CD4SPs 
was clearly faster and more effective in BDC2.5/B6 g7 
than in BDC2.5/NOD mice (Fig. 2 H and Fig. S3). 

The analysis of two additional TCR tg lines (OTII 
and OTI; Fig. 3), as well as previous publications de- 
scribing other TCR tg lines on the NOD background 
(3A9 and AI4; Lesage et al., 2002; Choisy-Rossi et al., 
2004; Liston et al., 2007; Serreze et al., 2008), revealed 
that these observations were neither peculiar to the 
BDC2.5 model nor limited to MHCII-restricted TCR 



Figure 3. Impaired positive selection of MHC-II and 
MHC-I restricted TCRs on the NOD background. Thymocytes 
were isolated from mice expressing the OTII or the OTI tg TCRs 
on B6 or (B6xNOD.H2 b ) Fl background (Fl.H2 b ). Shown are the 
distribution of thymocyte subsets (A), expression of CD4 and V(J 5 
tg TCR chain on CD4SP thymocytes from OTII mice (B), and 
number of total thymocytes (C). Horizontal lines indicate the 
mean. Data are representative of two independent experiments 
(n > 3 per genotype per experiment). 
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Figure 4. Higher DN to DP transition and saturation of selection niches in BDC2. 5/NOD mice. Single and mixed BM chimera experiments in 
(B69 7 xN0D) F1 hosts (F1.H29 7 ;A-H) and BrdU pulse-labeling experiment in BDC2.5TCR tg mice on B69 7 or NOD background (I). (A-F) Thymocytes were analyzed 
in Fl hosts, 8 wk after reconstitution with single (A-C) or mixed (D-F) BM cells isolated from BDC2.5 mice on B6 or NOD backgrounds. Shown are the distribution 
of thymocyte subsets (A and D), expression of the clonotypic BDC2.5TCR on CD4SP thymocytes and CD4 + splenocytes (B and E), and of CD3e on DP thymocytes 
(Cand F, thin lines indicate control staining). Mean thymic cellularities of single BM chimera: 53 ± 15, n = 8 (BDC2.5/B69 7 in F1.H29 7 ); 111 ±31, n = 8 (BDC2.5/ 
NOD in F1 .H29 7 ). (G) Expression of CD45.1 (NOD) and CD45.2 (B69 7 ) on DN, DP, and CD4SP thymocytes from Fl mixed chimeras. The ratios of BDC2.5/N0D:BDC2.5/B6 
donors are 1 :1 (top) and 1 :10 (bottom). (H) N0D/B6 ratios in thymic subsets of mixed chimeras, over a range of BDC2.5/N0D:BDC2.5/B6 donor ratios. Each line 
represents one individual chimeric mouse. (I) More efficient DN to DP transition in BDC2.5/N0D mice in noncompetitive settings. Number of BrdU + DP thymocytes 
generated per 1,000 BrdU + DN at indicated time points after a single BrdU injection. Error bars indicate standard deviation of the mean. Data are representative of 
two (A-C; n = 8 per type of chimera) and nine (D-H; n = 38 mixed chimera) and means of eight (I; n > 25 per genotype) independent experiments. 



tg models. To avoid backcrossing steps, we took advantage of 
the fact that the phenotypes observed in BDC2. 5/NOD mice 
were dominant in BDC2.5/F1(B6 b7 xNOD) mice (unpublished 
data) and compared OTII and OTI TCR tg mice on B6 versus 
(B6xNOD.H2 b ).Fl backgrounds. Fl mice expressing the 
tg OTII and OTI TCRs, as well as NOD mice expressing 
the 3A9 (in absence of the transgene-encoded HEL TCR 
ligand) and AI4 tg TCRs, displayed marks of impaired positive 
selection similar to those seen in BDC2. 5/NOD mice, including 
decreased SP/DP ratios (Fig. 3 A; Lesage et al., 2002, Fig. 6 a; 
Serreze et al., 2008, Fig. 3 A), higher frequencies of T cells 
expressing endogenous TCR genes (Fig. 3 B; Lesage et al., 2002, 
Fig. 3 c; Choisy-Rossi et al., 2004, Fig. 6 A), decreased TCR 
levels on DP thymocytes (Liston et al., 2007, Fig. 5 a; Serreze 
et al., 2008, Fig. 4), and higher total thymic cellularities (Fig. 3 C; 
Lesage et al., 2002, Fig. 6 c; Choisy-Rossi et al., 2004, Fig. 1 B). 

Greater DN to DP transition and saturation 

of positive selection niches in BDC2. 5/NOD mice 

These several observations indicated that variation in the efficacy 
of positive selection was responsible for some of the phenotypic 
differences shown in Fig. 2 A. Positive selection is a saturable 
phenomenon whose efficiency is dependent on the availability 

JEM Vol. 209, No. 3 



of TCR-specific selection niches (Huesmann et al., 1991). To 
determine whether T cell— intrinsic and/or niche effects were 
responsible, we set up single and mixed BM chimera experi- 
ments. When chimeras were generated by parallel transfer of 
BDC2.5/B6B 7 and BDC2.5/NOD BM into Fl hosts, chimeras 
reproduced the differences in thymic CD4/CD8 profiles, fre- 
quencies of BDC low CD4SP cells, and CD3 levels on DPs, 
which were previously seen in parental tg thymi (Fig. 4, A— C). 

In contrast, when BDC2.5/B6S 7 and BDC2. 5/NOD T cells 
matured together in the same thymic environment of mixed 
BM chimeras, and thus faced the same competition for positive 
selection niches, these differences were neutralized (Fig. 4, D— F). 
The proportion of DPs was restored to normal in BDC2.5/ 
B6 g7 -derived cells, and BDC2.5/NOD-derived cells matured 
effectively into CD4SPs, if anything with slightly higher effi- 
ciency. Accordingly, the frequency of BDC low cells and TCR 
levels on DP thymocytes were also normalized in mixed chi- 
meras (Fig. 4, E and F), therefore establishing that the positive 
selection defect in BDC2. 5/NOD mice was driven by a de- 
creased availability of positive selection niches. 

A second observation, which largely explains the difference 
in intraclonal competition between the parental tg strains, 
was made in these mixed BM chimera; we noticed that the 
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Figure 5. The greater DN to DP transition in adult and fetal thymi 
from BDC2.5/N0D mice is MHC independent. (A-C) Thymocytes were 
analyzed in indicated chimeras, 8 wk after reconstitution with mixed BM 
cells isolated from BDC2.5 mice on B6 and N0D.H2 b backgrounds. 
Expression of CD45.1 (NOD) and CD45.2 (B6) on DN and DP thymocytes, 
for a ratio of BDC2.5/N0D:BDC2.5/B6 donors of 1 :1 (A and B) and N0D/B6 
ratios in DN and DP subsets, over a range of BDC2.5/N0D:BDC2.5/B6 
donor ratios (C). (D-E) BDC2.5 tg FTOCs expressing H2 g7 or nonselecting 
H2 b MHC haplotypes were analyzed after 7 d of culture. Shown are the 
expression of CD4 and CD8 (D) and number of DP thymocytes in FTOCs 
of indicated genotypes (E). Horizontal lines indicate the mean. Data are 
representative of three (A-C; n > 5 per type of chimera) and two 
(D-E; n > 8 per genotype) independent experiments. 



relative proportion of thymocytes derived from BDC2.5/ 
NOD and BDC2.5/B6 g7 in the reconstituted thymic subsets 
was very different. At a 1:1 ratio of BM cells, chimeras 
showed the expected 1:1 ratio in the thymic DN population 
(Fig. 4 G, top). However, the representation of BDC2.5/ 
B6 g7 progeny was greatly reduced in the DP compartment, 
accounting for <5% of the whole population. As a result, 
CD4SP and peripheral T cell compartments were largely 
dominated by BDC2.5/NOD-derived cells and 10:1 ratios 
were required in the BM inoculum to get a balanced repre- 
sentation of BDC2.5/B6B 7 and BDC2.5/NOD cells at the 
DP stage (Fig. 4 G, bottom). On average, for various input 
ratios, the DN to DP transition was 10— 20x more efficient 
on the NOD background (Fig. 4 H). BrdU pulse-labeling 
experiments showed that these results held true in noncom- 
petitive settings, since much higher ratios of BrdU + DP ver- 
sus BrdU + DN cells were observed in BDC2.5/NOD thymi 
than in B6 g7 counterparts (Fig. 4 I). In hindsight, these results 
explained the greater number of DP thymocytes observed in 
adult (Fig. 2 B and Fig. 3 C; Lesage et al., 2002, Fig. 6 c; 
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Figure 6. In vivo thymic negative selection is functional in NOD mice. 

Mixed BM chimera experiments in (B69 7 xN0D)F1 hosts (BDCmi") or F1 hosts 
expressing a BDC mimotope at low (BDCmi L0W ) or high (BDCmi HIGH ) levels in 
thymic epithelial cells. (A) Thymocytes were analyzed in indicated chimeras, 
8 wk after reconstitution with mixed BM cells isolated from BDC2.5/B69 7 , 
BDC2.5/N0D, and polyclonal Fl mice (ratio 4:1:0.3, respectively). Expression of 
CD45.1 (NOD) and the clonotypic BDC2.5 TCR on thymic subsets. Numbers 
indicate the percentage of remaining BDC2.5/B6 (green) and BDC2.5/NOD 
cells (red) relative to the number of cells observed in BDCmr hosts. (B) Quan- 
tification of data shown in A. Percentage of remaining thymocytes relative to 
the number of cells observed in BDCmr hosts. Each dot represents one indi- 
vidual chimeric mouse. Horizontal lines indicate the mean. Data are represen- 
tative of three independent experiments (n > 8 per type of chimera). 



Choisy-Rossi et al., 2004, Fig. 1 B) and fetal (Zucchelli et al., 
2005) TCR tg thymi on the NOD background. 

This increased DN to DP transition in BDC2.5/NOD mice 
was not dependent on cognate MHC recognition, and thus not 
a consequence of differential negative selection, as BDC2.5/ 
NOD.H2 b -derived thymocytes also over-competed BDC2.5/ 
B6 counterparts at the DN to DP transition in mixed chimera 
settings in (B6xNOD.H2 b ).Fl hosts expressing the nonselecting 
H2 b MHC haplotype (Fig. 5, A and C), as well as in MHC- 
deficient hosts (Fig. 5, B and C). Symmetrically, at the fetal 
stage, the differential representation of DP thymocytes that was 
seen in BDC2.5/B6« 7 and BDC2.5/NOD FTOCs (Zucchelli 
et al., 2005) was also observed on the nonselecting H2 b MHC 
(Fig. 5, D and E), therefore excluding any contribution of nega- 
tive selection in this setting as well. 

In conclusion, these results pointed to the following scenario: 
a thymocyte-autonomous variation primarily affected the 
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Figure 7. B6 and NOD a(i-TCRs pro- 
mote opposite 78- and «|S-lineage fates. 

DN thymocytes (A-E) or splenocytes (F) were 
isolated from non-tg or BDC2.5 TCR tg mice. 
(A) Expression of CD3e, CD71, and CD5 on DN 
subsets (defined as shown in Fig. S4 A). (B) BrdU 
pulse-labeling experiment. Density plots show 
intracellular levels of incorporated BrdU in DN 
thymocytes at different times after a single BrdU 
injection. (C) Fold change/fold change plot de- 
picting the BDC2.5 tg/non-tg ratio of expression 
in B69 7 (x axis) and NOD (y axis) DN4 thymocytes. 
Blue lines, fold change of 2 (with the number 
of genes in the corresponding areas). (D) Heat 
map of 42 cell cycle-related transcripts down- 
modulated in BDC2.5/B69 7 DN4 thymocytes 
(relative expression is color-coded). Cycling (blast) 
and noncycling (small) DP and yh thymocytes 
from ImmGen datasets are included for com- 
parison. (E) Fold change versus mean expression 
plot depicting B69 7 versus NOD DN4 thymocyte 
subsets from BDC2.5 TCR tg (top) and polyclonal 
mice (bottom). Blue lines, fold change of 1.5. The 
number of 78 and a(3 signature genes above 
and below this threshold is indicated in green 
and red, respectively. (F) Proportion of 78-like 
CD4~CD8~ DNs on gated TCR-E5+ splenocytes 
from mice of indicated genotypes. Mean fre- 
quencies: 1 5.4 ± 4.9 (BDC2.5/B63 7 ) versus 0.6 ± 
0.2 (BDC2.5/NOD); 32.1 ± 1 1.3 (BDC2.5/B6) 
versus 3.8 ± 1.5 (BDC2.5/N0D.H2 b ). Data are 
representative of three (A; n > 9 per genotype), 
eight (B; n > 25 per genotype), two (C-E; n = 2 
per genotype per experiment), and three (F; n> 10 
per genotype) independent experiments. 



efficiency of the DN to DP transition, resulting in a greater num- 
ber of DP thymocytes in TCR tg thymi on the NOD back- 
ground, and thus more intraclonal competition for positive 
selection niches, ultimately leading to sluggish positive selection. 



In vivo negative selection 
is functional in NOD mice 

Because the previous analyses identi- 
fied the differential efficiency of the 
DN to DP transition as the primary 
cause of the phenotypes observed 
in TCR tg mice on the NOD back- 
ground, we reasoned that this same 
phenomena had likely confused the 
interpretation of the same models in 
conditions of in vitro and in vivo neg- 
ative selection (Lesage et al., 2002; 
Choisy-Rossi et al., 2004; Zucchelli 
et al, 2005; Liston et al, 2007; Serreze 
et al., 2008), putting in jeopardy the ar- 
guments that had been used to support a 
negative selection defect in NOD mice. 
18 **■ Therefore, we revisited the efficiency 

of negative selection in vivo on the 
NOD background, using mixed BM chimera in which niche 
effects were neutralized. We generated mixed BM chimeras 
into hosts that expressed a BDC mimotope at low (BDCmi LOW ) 
or high (BDCmi HIGH ) levels in thymic epithelial cells, using a 
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Figure 8. Selective attenuation of Erk1/2 phosphorylation downstream of NOD ap-TCRs, rather than global dampening of a|}-TCR signaling 
responses. Thymocytes (A-F) or splenocytes (G) were isolated from non-tg B69 7 and NOD mice or from mice expressing the BDC2.5 tg TCR on a fiog 7-deficient, 
nonselecting B6 and N0D.H2 b backgrounds, mixed in a 1:1 ratio and assessed forTCR-induced signal transduction in the same tube (CD45 congenic markers 
were used to distinguish B6 and NOD cells after stimulation), except in B where stimulation was performed in separate tubes. (A) Calcium flux induced in poly- 
clonal DP thymocytes by biotinylated primary antibodies (first arrow) and streptavidin cross-linking (second arrow). Third arrow, ionomycin control. Lower panels 
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tetracycline-regulated expression system (see Materials and 
methods). BDC2.5/B6 and BDC2.5/NOD BM cells were 
mixed (polyclonal Fl BM cells were also added to the mix of 
cells to avoid a massive clonal deletion in the reconstituted 
thymi) and used to reconstitute Fl control hosts or Fl hosts 
expressing the BDC2.5 mimotope. As expected, the neo- 
expression of the BDC2.5 mimotope induced the deletion of 
a significant number of BDC + thymocytes. On average, the 
number of BDC+ CD4SP was decreased by 55 ± 16% (n = 12) 
in hosts expressing the mimotope at low dose, whereas the 
level of deletion further increased to 95 + 6% (n = 10) for 
high doses of mimotope expression. Deletion took place 
mostly at the DP stage (Fig. 6) because both CD4SP and DP 
subsets were similarly affected by the deletion, whereas DN 
thymocytes remained mostly unaffected. Importantly, the rel- 
ative decreases in cell numbers observed in BDC2.5/B6 g7 and 
BDC2.5/NOD-derived thymocytes were similar at either 
dose of mimotope expression (Fig. 6). Thus, negative selec- 
tion, in vivo, is functional on the NOD background. 

Tg ap-TCRs promote opposite 78- and ap-lineage fates 
in B6 and NOD thymocytes 

Because the DN to DP transition is driven by the ctfi-TCR 
inTCR tg thymi (Trop et al., 2000), understanding the cellu- 
lar and molecular underpinnings of its differential efficacy on 
B6 and NOD backgrounds would potentially reveal how 
NOD genetic variation impinges on a(3-TCR signaling. 
Therefore, we analyzed the precursor DN pool (Fig. S4 A). 
TCRs were expressed on the cell surface at the transition be- 
tween the DN2a and DN2b stages, and at the same levels on 
the two backgrounds (Fig. 7 A). Expression of the TCR was 
associated with the up-regulation of CD5 and the prolifera- 
tion-associated marker CD71 (transferrin receptor). Whereas 
CD5 up-regulation was stronger on BDC2.5/B6 g7 DN thy- 
mocytes, likely indicative of stronger TCR signals, the ex- 
pression of CD71 was transient on these cells and largely 
down-modulated by the time cells reached the DN4 stage. In 
contrast, a higher degree of cell cycling was observed in the 
DN4 pool of BDC2.5/NOD mice, as assessed by the persis- 
tence of CD71 at this stage (Fig. 7 A) and the higher fraction 
of DN cells incorporating BrdU in pulse-labeling experi- 
ments (Fig. 7 B). Interestingly, although BDC2.5/NOD cells 
showed a rapid decrease in BrdU levels during the chase 



period, presumably as a result of their differentiation into 
DP cells, levels remained high in BDC2.5/B6 g7 DN cells, 
indicating that a significant fraction of them accumulated at 
the DN stage (Fig. 7 B). 

To characterize the transcriptional impact of the early ex- 
pression of the tg a(3-TCR at the DN stage on the two back- 
grounds, we then generated gene expression profiles from 
BDC2.5 tg and polyclonal DN2-3 and DN4 subsets. According 
to a fold change/fold change plot displaying the ratio of expres- 
sion in BDC2.5 tg versus non-tg DN4 thymocytes from B6 g7 
versus NOD mice (Fig. 7 C), the BDC2.5 TCR induced exten- 
sive gene expression changes in DN4 thymocytes. Among the 
genes repressed by the BDC2.5 TCR in DN4 thymocytes from 
B6 g7 mice, a Gene Ontology enrichment analysis identified a 
cluster of 42 cell cycle— related transcripts (Fig. 7 D). These genes 
were not repressed or were repressed to a much lower extent in 
BDC2.5/NOD DN4 thymocytes, consistent with the differen- 
tial degree of proliferation seen in the previous paragraph. 

To determine whether these broad transcriptional changes 
denoted a particular lineage, we used the most induced and re- 
pressed genes (603 gene probes, fold change > 3) to define the 
footprint of the BDC2.5 TCR on DN progenitors, and rank 
the cell populations analyzed by the Immunological Genome 
(ImmGen) Consortium (www.immgen.org; 267 populations 
covering lymphoid and myeloid lineages) according to the 
level of expression of these genes (Fig. S4 B). Strikingly, 16 of 
the 20 top ranked populations belonged to the •yS-lineage. 
This assignment was confirmed by highlighting a(3- and 
"y8-lineage commitment signatures (defined independently 
as described in Fig. S4 C) in Fig. 7 C. Most of the up-regulated 
genes belonged to the 78-lineage signature, whereas many of 
the down-modulated ones belonged to the a(3 commitment 
signature. These changes were consistent with earlier studies 
proposing that premature expression of tg a|3-TCRs in DN 
thymocytes commits a significant fraction of these progenitors 
to the 78-lineage (Bruno et al., 1996; Terrence et al., 2000). 
However, the off-diagonal placement, as well as the total 
number of gene probes up-regulated or down-regulated by 
the a(3-TCR on the two backgrounds (at an arbitrary thresh- 
old of twofold change, 1,933 versus 333 on B6 g7 and NOD 
backgrounds, respectively; Fig. 7 C), indicated that these 
transcriptional changes were preferentially induced on the 
B6 B? background, and were much more muted on the NOD 



display overlaid mean ratios. (B) Purified polyclonal DP thymocytes were left untreated (x) or stimulated (anti-CD3e/CD4) before isolation of whole-cell lysates 
and immunoblotting with antibodies specific for phosphotyrosine (pY) and Erkl/2. Arrowheads indicate presumed proteins based on their expected molecular 
size. (C) Erkl/2 phosphorylation induced in DP thymocytes from polyclonal or BDC2.5 TCR tg mice by anti-CD3e/CD28 antibodies or peptide-MHC tetramers. 
Cells were fixed at the times shown at left, after streptavidin or BDC2.5 tetramer addition (no cross-linking controls in the first row were 2 min after addition of 
primary antibodies alone). Last row, PMA positive control. Numbers indicate percentage of pErkl /2 + cells. (D) Representation of data shown in C as the fre- 
quency of p-Erkl /2 + cells against the mean fluorescence intensity (MFI) of p-Erkl /2 + cells. Time points (in min) are indicated in the circles, x, levels observed in 
absence of cross-linking. (E) Quantification of data shown in C and D, showing the frequency of pErk + NOD thymocytes at the peak of the response (t = 2 min) 
relative to the frequency of B6 responders. Horizontal lines indicate the mean. (F and G) Impaired Erkl/2 phosphorylation at multiple stages of T cell differentia- 
tion on the NOD background. Erk1/2 phosphorylation induced in CD25+ DN thymocytes from BDC2.5 TCR tg mice (F) and in polyclonal CD4+ and CD8+ spleno- 
cytes (G) by anti-CD3s/CD28 antibodies, in the same stimulation conditions as in C. Last row, PMA positive control. Numbers indicate percentage of pErkl/2 + 
cells. Data are representative of three (A; n = 3 per genotype), two (B; n = 2 per genotype), five (C-E; n = 17 per genotype for polyclonal stimulations and n > 5 
per genotype for BDC2.5TCR tg stimulations), three (F; n = 5 per genotype), and two (G; n = 4 per genotype) independent experiments. 
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background. This differential magnitude reflected intrinsic 
properties of B6 and NOD a(3-TCR signaling networks and 
not pre-TCR or generic strain-specific differences because 
they were not observed in polyclonal mice (Fig. 7 E). Accord- 
ingly, B6 but not (or to a much lower extent) NOD a(3-TCR 
signalosomes induced most of the transcripts found in thymic 
78T cells, which include the 78-lineage transcription factor 
Soxl3 (Ciofani and Zuhiga-Pfliicker, 2010), Slprl, which 
promotes cellular egress from the thymus, and genes that are 
strictly expressed in the 78 branch of thymocytes, like Etv5, 
Blk, Cpa3, AM, Rbpms, and Tox2 (Fig. 7 C and Fig. S4 D). 

To determine whether these commitment programs were 
functional and effectively diverted DN progenitors toward 
the 78-lineage, we directly compared the transcription factors 
induced by the BDC2.5 TCR in DN4 thymocytes to the 
ones equipping polyclonal 78T cells (Fig. S4 E, top). The 
BDC2.5 TCR induced the expression of >70% of the tran- 
scription factors expressed in 78T cells. To identify the 
molecular events initiating 78T cell fate decision downstream 
of B6 a|3-TCR signalosomes, we analyzed the DN2-3 thy- 
mocyte subsets that corresponded to the earliest stages of 
expression of the a(3-TCR (Fig. S4 E, bottom). Of 188 78 
transcription factors, only 21 were induced by the BDC2.5 
TCR, including previously unknown 78 fate regulators like 
Etv5, Tox2, Rbpms, and Sox5, as well as the already identified 
Soxl3 and Id3 factors. In accordance with the extensive array of 
78 transcription factors induced downstream of B6 a(3-TCR 
signalosomes, CD4~CD8~TCR|3 + peripheral cells, described 
earlier as 78-like T cells in TCR tg models (Bruno et al., 
1996; Terrence et al., 2000), constituted a substantial fraction 
of the T cell pool in BDC2.5/B6« 7 mice, and were 10-20X 
more abundant than in NOD mice, both on selecting and 
nonselecting backgrounds (Fig. 7 F). 

Therefore, in addition to promoting the expression of 
78 fate-associated genes, B6 but not NOD a|3-TCR sig- 
nalosomes induced a full 78 maturation program, associated 
with thymic export and colonization of peripheral lym- 
phoid organs. In contrast, NOD DNs expressed more of the 
a|3-lineage genes and the proliferative program character- 
istic of ot(3-committed DN thymocytes. This differential 
branching likely explains the much higher efficacy of the DN 
to DP transition in BDC2.5/NOD mice. 

Selective Erk1/2 defect downstream of ap-TCRs in NOD T cells 

What is the biochemical nature of the signaling differences 
driving these contrasting cell fates downstream of B6 and 
NOD a(3-TCR signalosomes at the DN stage? As a(3-TCR 
rather than pre-TCR signaling was affected (Fig. 7 E), we an- 
ticipated that the same signaling defects would be present at 
other stages where the a|3-TCR is physiologically expressed 
in normal mice, in particular at the DP stage. In addition, be- 
cause current models of ap/78 commitment propose that 
the differential strength of signals downstream of pre-TCRs 
and 78-TCRs drive these opposite commitments, we quan- 
tified signal transduction downstream of B6 and NOD 
a(3-TCRs along several known pathways in DP thymocytes. 
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We first analyzed calcium mobilization, which is depen- 
dent on a functional proximal signaling module (Lck— Zap70— 
Lat— PLC7-I). Intracellular calcium concentrations reached 
very similar levels in polyclonal B6 B? and NOD DP thymo- 
cytes after TCR stimulation (Fig. 8 A). We then quantified 
global tyrosine phosphorylation events by immunoblotting 
extracts from TCR-stimulated DP thymocytes. No obvious 
defects in the patterns of phosphorylated proteins were de- 
tected in NOD cells (Fig. 8 B). These two results excluded the 
hypothesis of a general dampening of NOD TCR signaling 
pathways as the cause of the opposite ap/78 commitments 
observed in TCR tg mice on NOD and B6 backgrounds. We 
also analyzed the activation of the Erkl/2 MAPK, one of the 
key pathways known to be differentially required during 
TCR-dependent checkpoints in thymocyte differentiation 
(Ciofani and Zuniga-Pflucker, 2010; Teixeiro and Daniels, 
2010). The fraction of cells showing Erk phosphorylation 
(Fig. 8, C and E), as well as the levels of phosphorylation 
(Fig. 8 D), were much lower on the NOD than on B6 g7 back- 
ground for all time-points tested after either anti-TCR antibody 
stimulation or more physiological antigen-specific stimulation 
by peptide-MHC tetramers. 

Of note, this defect was not only present in DP thymocytes 
but also at the other stages where the a(3-TCR was expressed, 
including in anti-TCR— stimulated BDC + DN2/3 thymocytes 
(Fig. 8 F), which corresponded to the earliest stages of expres- 
sion of the a(3-TCR in BDC2.5 mice (Fig. 7 A), as well as in 
anti-TCR— stimulated peripheral CD4 + and CD8 + T cells from 
non-tg mice (Fig. 8 G) . Together, these results pointed to a se- 
lective defect in the Erkl/2 signaling module downstream of 
a(3-TCRs on the NOD background. 

DISCUSSION 

Using several TCR tg mouse models to study thymocyte se- 
lection on the NOD background, we show that a marked 
blockade of positive selection of clonotypic T cells occurs in 
thymi from TCR tg mice on the NOD background, a block- 
ade caused by higher numbers of DP thymocytes competing 
for positive selection niches. These higher numbers themselves 
reflect an enhanced DN to DP transition, which in turn results 
from the quasi-absence of deviation to the 78-hneage, usually 
imparted by early expression of tg TCRs (Fig. 9). These obser- 
vations held true on nonselecting backgrounds, a result which, 
by excluding any contribution of MHC-dependent mecha- 
nisms, formally invalidates our earlier model according to 
which increased DP thymocyte numbers in BDC2.5/NOD 
result from impaired clonal deletion mechanisms (Zucchelli 
et al., 2005; Holler et al., 2007). Accordingly, using an in vitro 
clonal deletion model based on isolated, preselected, and 
monoclonal DP thymocytes on one hand and in vivo negative 
selection conditions in mixed BM chimera on the other hand, 
we found that ct(3-TCR complexes equipped with either B6 
or NOD signaling networks show very similar signaling po- 
tency and commit DP thymocytes to apoptosis with the same 
efficiency upon encounter of negatively selecting ligands. Thy- 
mic negative selection is therefore functional in NOD mice. 
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Figure 9. Proposed model, outlining the differential signaling abilities of NOD and B6 ap-TCRs to promote ctp- and 78-lineage commit- 
ment, when expressed at the DN stage. As a result of the high levels of tg af3-TCR expressed on DN thymocytes in TCR tg mice, a|3-TCRs mimic 
78-TCRs and divert a significant number of DN progenitors toward the 78-lineage on the B6 background (left). Because of a selective defect in the Erk1/2 
signaling module, the same tg TCRs wired to NOD signalosomes are perceived like pre-TCRs, rather than 78-TCRs, leading to amplification and abundant 
differentiation into DPs. As a consequence, thymi from NOD TCR tg mice show higher numbers of DP thymocytes than their B6 counterparts, and relatively 
inefficient positive selection, promoting compensatory rearrangements of endogenous TCR-a loci which rescue some measure of positive selection. 
This differential decision at the ap/78 commitment checkpoint explains most of the phenotypes reported in TCR tg NOD mice, among which the impaired 
selection of TCR tg thymocytes as a result of the saturation of positive selection niches, the higher frequency of receptor editing and clonotype' 0 ™ T cells, 
and the decreased TCR levels on DP thymocytes. Of note, mostT reg cells emerge from the clonotype' 0 ™ fraction in TCR tg models and the higher frequency 
of clonotype' 0 ™ T cells in BDC2.5/N0D mice results in higher T reg/T conv cell ratios, likely explaining decreased diabetes incidence on this backgrounds. 



As purposely emphasized in the Results section, the hall- 
marks of the divergence at the afJ/78 commitment check- 
point seen in BDC2.5/NOD mice are similarly present in 
the other TCR tg models for which negative selection de- 
fects were reported on the NOD background. It thus seems 
very likely that, as observed with the BDC2.5 model (and 
also extended here to the OTII and OTI models), an in- 
creased DN to DP transition also leads to the saturation of 
selection niches in the 3A9 and AI4 systems (indeed, differ- 
ential DN to DP transition was seen in the mixed BM chi- 
meras of Lesage et al., 2002). Together with a differential 
efficiency of 78 deviation (and indeed, much lower propor- 
tions of 78-like TCR-|3 + DN T cells were recently reported 
in 3A9 mice on the NOD compared with BIO and BALB/c 
backgrounds; Dugas et al., 2010), this phenomenon would 
have obscured an efficient thymic negative selection. Although 
these predictions await verification, interpretation of the 
published data in light of our current results strongly suggests 
that the same mechanisms are at play in these models. 

Because of the early TCR expression in TCR tg models 
and the higher affinity of TCR-a for TCR-(3 compared with 
pTa (Trop et al., 2000), a(3-TCR rather than pre-TCR sig- 
nalosomes drive the DN to DP transition in TCR tg models. In 
accordance with earlier hypotheses that tg a(3-TCRs divert a 
significant fraction of DN thymocytes toward the 78-lineage 
(Bruno et al., 1996; Terrence et al., 2000; Egawa et al., 2008), 
transcriptome analyses formally demonstrated that B6 a(3-TCR 
signalosomes indeed promote cell-autonomous 78 commitment 



programs of thymic progenitors, early maturation, and thymic 
egress before the DP stage. Therefore, these cells do not con- 
tribute to the DP compartment, explaining the lower DP cel- 
lularities and the efficient selection of the BDC2.5 clone seen 
on the BDC2.5/B6 g7 background. Of note, this 78 diversion in 
TCR tg models was reported on several genetic backgrounds 
(including B6, B10, and BALB/c; Terrence et al, 2000; Dugas 
et al., 2010). Therefore, it does not reflect a peculiarity of the 
B6 background. In contrast, the signaling properties of NOD 
a(3-TCR signalosomes stand out as they resemble those of 
pre-TCR complexes by promoting a(3 rather than 78 com- 
mitment of DN thymocytes (Fig. 9). Contrary to most other 
TCR tg models on various genetic backgrounds, TCR tg mice 
on the NOD background therefore show an unusually normal 
thymic differentiation, with cell numbers and subset frequen- 
cies that are much closer to a non-tg thymus. 

How to explain then that NOD and B6 ct(3-TCR signalo- 
somes promote equivalent outcomes at the DP stage (cell death 
by negative selection) but opposite cell fates (a(3 versus 78 com- 
mitment) at the DN stage? Based on the quantitative models of 
afJ/78 commitment, 78 commitment requires stronger TCR 
signals than a(3 commitment and defects in TCR signaling or 
ligand availability were shown to convert 78 into a (3 commit- 
ment signals (Ciofani and Zuniga-Pfliicker, 2010). However, 
because quantitative measurements of TCR signaling did not 
reveal a global dampening or proximal defects in NOD TCR 
signaling pathways, a generic decrease in TCR signaling strength 
is not likely to account for the differential afi/78 commitment 
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decisions observed in TCR tg models on NOD and B6 back- 
grounds. In contrast, the selective defect in Erkl/2 activation on 
the NOD background makes a particularly attractive candidate, 
as Erkl/2 and the Id3 levels they control were shown to be cru- 
cial regulators of 01(3/78 commitment, with the adoption of the 
■78 fate being associated with greater Erkl/2 phosphorylation 
(Hayes and Love, 2002; Haks et al., 2005; Hayes et al, 2005). 
Because tg ct(3-TCRs are expressed at high levels on DN thy- 
mocytes compared with pre-TCR complexes in polyclonal 
mice, we propose that such high receptor density conveys 
strong, 78-like Erkl/2 signals that instruct 78-lineage differen- 
tiation on the B6 background. In contrast, on the NOD back- 
ground, high levels of tg a(3-TCRs do not elicit such high 
Erkl/2 signals and behave like pre-TCRs. Of note, a recent 
study showed that pre-TCR— mediated a|3 commitment is less 
efficient in NOD than in B6 polyclonal mice (Feng et al., 201 1). 
Because Erkl/2 defects are likely present downstream of NOD 
pre-TCRs and because pre-TCRs are expressed in limiting 
amounts on the cell surface of DN thymocytes, contrary to the 
highly expressed a|3-TCRs in TCR tg settings, we anticipate 
that the same defect might hamper (rather than promote) the 
pre-TCR— driven DN to DP transition in polyclonal settings, a 
prediction which agrees with the aforementioned study (Feng 
et al., 2011). Therefore, the results observed here in TCR tg 
settings are not direcdy transposable to polyclonal settings, as a 
result of the very different density of receptors driving the DN 
to DP transition in TCR tg versus polyclonal settings. 

Erkl/2 kinases were shown to be dispensable during thy- 
mocyte negative selection (McGargill et al., 2009). Therefore, 
this selective defect fits well with the observation that clonal 
deletion and associated transcriptional responses are intact 
on the NOD background. Similarly illustrating that not all 
functional outcomes are equally dependent on the Erkl/2 
module, the conditional deletion of Erk2 in peripheral CD8 + 
T cells recently revealed its selective requirement for T cell 
survival but not for T cell activation or effector functions 
(D'Souza et al., 2008). The depressed phosphorylation of 
Erkl/2 is also present in mature NODT cells. The functional 
relevancy of this defect in peripheral T cells forTID will 
require further investigation, but it is noteworthy that T cell 
defects in the phosphorylation of Erkl/2 were recently 
reported in two other autoimmune conditions, in lupus 
(Gorelik and Richardson, 2010) and rheumatoid arthritis 
(Singh et al.,2009). 

Erkl/2 kinases are required for positive selection of con- 
ventional T cells (T conv cells; Fischer et al., 2005) and agonist- 
selected T reg cells (Willoughby et al., 2007). Surprisingly, 
however, the reduced Erkl/2 signaling on the NOD back- 
ground has no discernable effect on the positive selection of 
thymocytes expressing the fixed BDC2.5 TCR in mixed chi- 
mera settings. If anything, selection is actually more efficient 
on the NOD background. One might speculate that, depend- 
ing on the affinity of each particular TCR for its positively 
selecting ligands, such a selective dampening of Erkl/2 re- 
sponses in NOD thymocytes might favor, impair, or not affect 
its selection by shifting the Erkl/2 signal inside or outside 
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the optimal window of selection. Thus, when considering a 
polyclonal repertoire, one might predict that the affinity 
threshold for positive selection, and thereby overall affinities 
of T conv TCRs for self-MHCp, may be higher in NOD 
mice. Because Erkl/2 kinases are dispensable for effector 
functions of peripheral T cells (D'Souza et al., 2008), this se- 
lective defect would not attenuate the higher reactivity of pe- 
ripheral T cells, in agreement with the enhanced proliferative 
T cell responses recently reported in NOD mice (D'Alise et al., 
2008). Conceivably, by narrowing the set of self-peptides that 
would reach the affinity and pERK thresholds required for 
T reg cell selection, this defect might also explain recent re- 
sults showing that, despite their apparent normal thymic de- 
velopment, TCR diversity within the NOD T reg cell 
compartment is markedly restricted (Ferreira et al., 2009). 
Supporting these hypotheses, the positive selection of invari- 
ant NKT (iNKT) cells, which is also Ras-Mapk-Egr depen- 
dent (Hu et al., 2011), is impaired in NOD mice that have 
few thymic iNKT cells (Gombert et al., 1996). Because iNKT 
cells express an invariant TCR, no repertoire compensation 
could counterbalance the Erkl/2 defect, potentially explain- 
ing their decreased numbers in NOD mice. 

In conclusion, the cellular and molecular dissection of 
thymic selection in NOD mice using multiple TCR tg mod- 
els yields two key conclusions. First, thymic negative selection 
does not appear to be at the root of autoreactivity in this 
model. Second, instead of a global dampening, NOD TCR 
signaling networks show a selective defect in the Erkl/2 sig- 
naling module. We therefore propose the working hypothesis 
that the emergence of autoreactivity might be favored in 
NOD mice by a shifted positive selection window and an 
overall higher reactivity of peripheral T cells for self-MHCp. 
Although new genetic tools will be required to directly test 
this model, it constitutes a new lead to the broader question of 
the mechanisms by which central tolerance failures lead to the 
emergence of autoreactive T cells in autoimmune diseases. 

MATERIALS AND METHODS 

Mice. All mice were maintained in specific pathogen-free facilities at Harvard 
Medical School (Institutional Animal Care and Use Committee, 02954). 
BDC2.5 TCR tg mice on C57BL/6.H2B 7 and NOD background were pre- 
viously described (Katz et al., 1993). BDC2.5/NOD.Rag _/_ and BDC2.5/ 
B6 s7 .Rag _/_ mice were derived by crossing the previous lines with NOD. 
Rag _/_ and B6 g7 .Rag _/_ lines, respectively, and diabetes occurrence was 
prevented by neonatal i.p. injection of total spleen cells from (B6 g7 xNOD) Fl 
mice, as described previously (Gonzalez et al., 2001). Similar breeding 
strategies were used to obtain BDC2.5 TCR tg lines on Rrt^-sufficient or - 
deficient B6 and NOD.H2 b background (NOD.H2 b line was previously 
described; Wicker et al., 1992). TCR-cr'-/NOD, MHC-deficient (IA0-'- 
P2m -/ ~), OTII, and OTI.Rag -7- mice were previously described (Philpott 
et al., 1992; Grusby et al., 1993; Barnden et al., 1998; Clarke et al., 2000). 
OTII/Fl.H2 b and OTI.Rag _/ ~/Fl.H2 b mice were derived by crossing the 
latter with NODH2 b and NODH2 b .Rag _/ ~ mice, respectively. For BM chi- 
meras, polyclonal and BDC2.5TCR tgmice on (B6s 7 xNOD) F1.H2S 7 back- 
ground were generated. NOD mice expressing the BDC mimotope were 
generated by crossing the tetO-BDCmi mouse line (described in the next 
section; on the NOD background) with CII-tTA mice (Witherden et al., 
2000; expressing the transcriptional activator under the control of the MHC 
Ea promoter, driving the expression in thymic epithelium; backcrossed on 
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the NOD background for >12 generations). Unless otherwise specified, all 
mice were analyzed between 4 and 6 wk of age. 

Generation of tetO-BDCmi mice. To generate this mouse line encoding 
a tetracycline (tet) -inducible BDC2.5 mimotope (BDCmi#l, sequence = 
RTRPLWVRME), the pTiB construct was generated from the pTIM vector, 
as previously described (van Santen et al., 2004). This vector consists of seven 
tetO sequences, a CMV minimal promoter, rabbit [3-globin intron 2, the 
mouse class II— associated invariant chain (Ii) peptide region in which amino 
acids 88—89 have been replaced by an MCC peptide, and the rabbit f}-globin 
3'-UTR with a poly(A) + RNA site. The pTiB plasmid was constructed via a 
two-step overlap PGR strategy replacing the MCC peptide with the BDC2.5 
mimotope peptide. A Pvul— Xhol fragment derived from pTiB was injected 
into NOD fertilized eggs and transgene positive founders were selected. 

Cell purification and sorting. DN thymocytes were purified by CD4/CD8 
complement-mediated depletion, as previously described (Laurent et al., 2004), 
and subdivided by cell sorting using an ARIA-II (BD). A dump channel was used 
to exclude CD19, CDllc, CDllb, CD49b, TCR-75, and GRl-positive cells. 
For polyclonal DN populations, CD3e was included in die dump channel. For 
clonal deletion assays and biochemistry, DP thymocytes were purified by mag- 
netic CD8 positive selection (Invitrogen; purity >95%). For microarrays, cells 
were double sorted (purity > 99.9%) directly into Trizol. 

Flow cytometry. All fluorochrome-conjugated and biotinylated antibodies 
were purchased from BD, eBioscience, or BioLegend, except for the biotinyl- 
ated antibody directed against the BDC2.5 TCR (Kanagawa et al., 2002), 
which was produced in the laboratory, and for biotinylated I-A g7 /BDC2.5 
mimotope monomers (AHHPIWARMDA), which were obtained from the 
National Institutes of Health Tetramer Facility. Cell events were collected using 
an LSR-II cytometer (BD) and analyzed with Flowjo software (Tree Star). 

Microarrays. RNA was prepared from double-sorted cell populations using 
Trizol, as previously described (Yamagata et al., 2004). RNA was amplified for 
two rounds (MessageAmp aRNA; Invitrogen), biotin labeled (BioArray High 
Yield RNA Transcription Labeling; Enzo Life Sciences), and purified with the 
RNeasy Mini kit (QIAGEN). The resulting cRNAs were hybridized to 
GeneCliip Mouse Genome M1.0 ST chips arrays (Affymetrix) . The image reads 
were processed through Affymetrix software to obtain raw .eel files, background 
corrected and nomialized using the RMA algoridnn implemented in the Expres- 
sion File Creator module from the GenePattem software package (Reich et al., 
2006). The cell populations analyzed were generated in duplicate and repli- 
cates were averaged. Data were analyzed with the Multiplot module from 
GenePattem. Microarray data are available from the National Center for Bio- 
technology Infomiation/GEO repositoiy under accession no. GSE34936. 

In vitro clonal deletion assay. Purified DP diymocytes from B6 and NOD 
background were mixed in a 1:1 ratio and co-cultured for various times with 
CD3/CD49b-depleted splenocytes from TCR-a _/ 7NOD mice (ratio DP/ 
APC of 1:2.5) that were previously pulsed for 2 h with graded concentrations of 
mimotope peptides (BDCmi#l: RTRPLWVRME Qudkowski et al., 2001]; 
#2: RTLALWVRME; #3: KVLAVWVRMME) in serum-free medium or 
were incubated in medium only (for the unstimulated condition). The frequency 
of remaining viable cells was measured using 0.1 j-ig/ml DAPI, Annexin V 
(according to the manufacturer's instructions; BD), and SSC/FSC plots. 

FTOCs. Fetal thymus lobes were dissected from E15.5 embryos. Organs 
were cultured on sponge-supported filter membranes (Gel Foam surgical 
sponges, Pharmacia; Supor450 membrane, 0.45-um pore size; Pall Gelman 
Laboratory, VWR) at an interphase between 5% C0 2 -humidified air and 
complete culture medium. Cell suspensions were prepared from FTOC cul- 
tures at day 7 by mechanical disruption of the lobes. 

BrdU labeling. 1.6 nig BrdU (Sigma-Aldrich) was injected i.p. at the indicated 
time points before harvest. Thymocytes were fixed, pemieabilized, and stained 
using the APC BrdU flow kit (BD) per the manufacturer's instructions. 



BM chimeras. Lethally irradiated host mice (1,000 rad) were injected i.v. 
with 5 x 10 6 CD3/CD5/CD49b-depleted BM cells isolated from the femurs 
of donor mice. Mixed chimeras were performed by mixing BDC2.5/B6 and 
BDC2.5/NOD donor cells at indicated ratios. For analysis of in vivo clonal 
deletion, NOD mice carrying the tetO-BDC2.5mi transgene (generated as 
described in Generation of tetO-BDCmi mice) and the CII-tTA transgene 
(Witherden et al., 2000) were crossed with B6 g7 mice to generate 
(B6 g7 xNOD) Fl mice carrying either one or the two transgenes. The com- 
bination of the two transgenes drove a very strong expression of the BDCmi 
that resulted in an almost complete deletion of BDC2.5 T cells. Two other 
settings providing milder expression of the mimotope were therefore estab- 
lished, one that gave very low expression of the mimotope (BDCmi LOW ; 
tetO-BDC2.5 construct in absence of the CII-TA tg, corresponding to a low, 
TA-independent, leaky expression of the tetO promoter; Lewandoski, 2001) 
and a second that gave higher expression of the mimotope (BDCmi HIGH ; 
tetO-BDC2.5 construct combined with CII-TA tg, in the presence of 20 ug/ml 
doxycycline [Sigma-Aldrich], provided in drinking water for the whole time 
of the reconstitution). Doxycycline-supplemented water was changed three 
times per week. Chimeras were analyzed 8 wk after reconstitution. 

Calcium flux measurement. B6 and NOD total thymocytes were resus- 
pended in a 1:1 ratio before loading with the calcium indicator Indo-1 (In- 
vitrogen) as previously described (Wang et al., 2005). Baseline Indo-1 
fluorescence was recorded at 37°C before the addition of 5 fig/ ml of bio- 
tinylated stimulatory antibodies, 10 u.g/ml streptavidin, and 1.4 uM iono- 
mycin. Calcium responses were analyzed with Flowjo kinetics platform. 

Biochemistry. 5.10 6 DP thymocytes were resuspended in a volume of 100 ul 
and stimulated for 2 min at 37°C with 5 u.g of biotinylated anti-CD3, 5 ug anti- 
CD4, and 7 ug avidin. After lysis in TNE buffer (50 mMTris, 1% Nonidet 
P-40, and 20mM EDTA) supplemented with protease and phosphatase inhibi- 
tors, lysates were subjected to immunoblot analysis. Antibodies against phos- 
photyrosine and Erkl-2 were purchased from Cell Signaling Technology. 

Phospho-flow analysis. B6 and NOD total thymocytes were resuspended 
in a 1:1 ratio in medium containing 5 ug/ml of the biotinylated stimulatory 
antibodies and incubated for 2 min at 37°C before addition of 20 ug/ml 
streptavidin. At various times after cross-linking (or addition of 10 ug/ml 
I-A g7 /BDC2.5 tetramer), cells were fixed, pemieabilized, and stained as pre- 
viously described (Adachi and Davis, 2011). PMA controls were fixed 3 min 
after addition of 10 ng/ml PMA. 

Online supplemental material. Fig. SI shows the flow cytometry isolation 
strategy of preselected B6 and NOD DP thymocytes and kinetic microarray 
analysis of thymic negative selection programs on the two backgrounds. Fig. S2 
shows the flow cytometry gating strategy used to identify clonotype low T cells in 
BDC2.5 TCR tg mice. Fig. S3 shows the flow cytometry gating strategy used to 
generate the BrdU plot shown in Fig. 2 H. Fig. S4 shows the flow cytometry 
gating strategy and microarray datasets used to isolate and characterize DN diy- 
mocyte subsets in BDC2.5 TCR tg mice. Online supplemental material is avail- 
able at htt^://www.jem.org/cgi/content/fiill/jem.20112593/DCl. 
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